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ABSTRACT 

We predict the rate at which Gamma-Ray Burst (GRB) afterglows should be 
^ ■ detected in supernova searches as a function of limiting flux. Although GRB afterglows 

are rarer than supernovae, they are detectable at greater distances because of their 
higher intrinsic luminosity. Assuming that GRBs trace the cosmic star formation 
history and that every GRB gives rise to a bright afterglow, we find that the 
average detection rate of supernovae and afterglows should be comparable at limiting 
magnitudes brighter than K = 18. The actual rate of afterglows is expected to be 
Q\ • somewhat lower since only a fraction of all 7-ray selected GRBs were observed to 

have associated afterglows. However, the rate could also be higher if the initial 7-ray 
^r) ■ emission from GRB sources is more beamed than their late afterglow emission. Hence, 

current and future supernova searches can place strong constraints on the afterglow 
! appearance fraction and the initial beaming angle of GRB sources. 

q . Subject headings: gamma rays: bursts 



1. Introduction 

Since their discovery in the late 1960's (Klebasadel et al. 1973) through early 1997, 
Gamma-Ray Bursts (GRBs) had defied all attempts to determine their distance scale conclusively. 
The Burst And Transient Source Experiment (BATSE) on board the Compton Gamma-Ray 
Observatory (GRO) showed that the burst population is highly isotropic (Meegan et al. 1993; 
Briggs et al. 1993), suggesting that bursts occur at cosmological distances or in an extended 
Galactic halo. Moreover, the cumulative number counts of faint bursts deviated from that of a 
uniform distribution of sources in Euclidean space and flattened at faint fluxes, consistent with the 
expected effect of a cosmological redshift (Fishman & Meegan 1995, and references therein). Last 
year, with the advent of the BeppoSAX satellite (Boella et al. 1997), it became possible to localize 
GRB sources to within an arcminute on a timescale of hours. Such fast, accurate localizations were 
quickly followed by the detection of delayed X-ray (Costa et al. 1997), optical (van Paradijs et al. 
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1997), and radio (Frail et al. 1997) counterparts to GRB sources. In particular, Fell and Mgll 
absorption lines were detected at a redshift z = 0.835 in the spectrum of the optical counterpart 
to GRB970508 (Metzger et al. 1997), demonstrating conclusively that this burst occurred at a 
cosmological distance with a redshift z > 0.835. The isotropy of the burst population and the 
flattening of their number counts, taken in combination with the fact that the first confirmed 
redshift for an optical counterpart is high, provides strong evidence that GRB sources are located 
at cosmological distances. 

Most plausible GRB models involve either the collapse of a single massive star (e.g. Usov 
1992; Woosley 1993; Paczyhski 1998), or the coalescence of two compact objects - two neutron 
stars or a neutron star and a black hole - in a binary system (e.g. Paczyhsky 1986; Eichler et 
al. 1989; Narayan et al. 1992; Mochkovitch et al. 1993; Rees 1997). Since the lifetime of these 
progenitors is short compared to the Hubble time at a redshift z 5, the cosmic GRB rate should 
simply be proportional to the star formation rate at these redshifts, without any appreciable delay 
due to the finite progenitor lifetime. The cosmic rate of massive star formation rate has been 
determined from the U and 5-band luminosity density in Hubble Deep Field (Madau et al. 1996; 
Madau 1996; Madau, Pozzetti, & Dickinson 1997; Madau 1997). The inferred star formation rate 
ps(z) can then be converted to a GRB explosion rate Rgkb(z), based on the requirement that 
the latter would fit the observed number count distribution of 7-ray selected GRBs (Wijers et al. 
1997). 

Cosmological GRBs are at least 10 times rarer than Type II supernovae (SNell) - possibly 
even ~ 10 6 times rarer if GRBs occur primarily at high redshifts following the cosmic star 
formation history (Wijers et al. 1997). However, at peak luminosity, the GRB afterglows are 
~ 10 3 -10 4 times brighter than SNell. In Euclidean space, this would imply that GRBs are 
detected from a volume bigger by a factor ~ (IO 4 ) 3 / 2 = 10 6 , roughly canceling out the factor by 
which they are rarer than supernovae. Hence we expect that at some relatively bright limiting 
flux, the rate of afterglow detections should become comparable to that of SN detections. Current 
and future supernova searches should provide information about the fraction of GRBs which 
produce detectable afterglows. The statistics of bursts in 1997 for which afterglows could have 
been identified implies that this fraction is of order tens of percent (e.g., Castro-Tirado 1998). On 
the other hand, there could also be a population of afterglows without a GRB precursor. This 
would occur if the source emits a jet from which the 7-ray emission is more beamed than the 
subsequent optical afterglow radiation due to the deceleration of the jet by the ambient gas and 
the corresponding decline in its relativistic beaming with time (Rhoads 1997). A jet geometry 
would imply a higher rate of afterglow detections in supernova searches. 

In this Letter, we predict the detection frequency of GRB afterglows as a function of limiting 
flux at various observed wavelengths, and compare this rate with the analogous predictions for SNe 
Type la and Type II at high redshifts. We assume throughout a flat, Q = 1, A = 0, cosmology, 
with a Hubble constant Hq = 50 km s _1 Mpc -1 . 
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2. Input Parameters 

The inferred energy release in cosmological GRBs and their afterglows is comparable to the 
binding energy of a neutron star, ~ 10 53 ergs. We therefore assume that GRBs, like supernova 
explosions, have evolved stars as their progenitors. Thus, the cosmic GRB and SN rates should 
both trace the cosmic star formation history. For progenitors with an effective lifetime r, the 
corresponding event rate at time t > r will be proportional to the star formation rate at time 
(t — r). Following Dwek (1997), we get for the rate of supernovae Type la (SNela) per comoving 
volume, 

u , I 3 Is 6 'H m b)dm h J^p s [t-T(pm b )]f(p)dp 

#SNIa(t) = ™; T7 \ t , I 1 ) 

J mq>{m)dm 

where p s is the stellar mass formed per unit comoving volume per unit time, m denotes masses in 
solar units, (f)(m^) is the initial mass function (IMF) for binary systems with total mass in the 
range rri\ < < m u , f(p) = 24/i 2 is the distribution of secondary-to-total-binary-mass ratios 
fi = m/nib, r(m) is the lifetime of a progenitor star of mass m, and (3 < 1 is the SNIa rate 
amplitude, an adjustable parameter. We find that (3 = 0.05 yields the best fit to the local and 
high-redshift supernova observations, and adopt this value. The lower limit in the second integral 
of the numerator is fi m = max[m(t)/mb, (m-b — 8)/mb], where m(t) is the mass of stars which are 
turning off the main sequence at time t, i.e. the inverse of r(m) (Greggio & Renzini 1983). This 
lower limit take into account the fact that stars more massive than the turnoff mass should not 
be included in the integral. We assume a Salpeter (1955) IMF, (p(m) oc m -2 ' 35 , with m\ = 0.1 and 
m u = 125. For Type II supernovae, the progenitor lifetime is negligible compared to the Hubble 
time, and so we take r = 0, 

p B (t)Uw <t>{rn)dm + (1 - g) / 8 16 c/>(m)dm] 
J mi fn4>\ m )dTn 

The second term in the numerator accounts for single stars and members of binary systems which 
do not lead to Type la supernovae. 

If GRBs have stellar progenitors, then their rate also traces the star formation rate. This 
scenario was investigated by Wijers et al. (1997), who also assumed that the progenitors are 
short-lived compared to the Hubble time (r = 0), and derived a best-fit constant of proportionality 
between the GRB rate -Rgrb and the star formation rate p s . Given our assumptions about the 
IMF, their result translates to 

R G RB(t) = 7.1 x 10- 7 #snii(£). (3) 

Note that this rate is substantially lower (by a factor ~ 150) than the best-fit rate derived by 
assuming a non-evolving burst population, i?GRB ~ 10~ 4 -Rsnii- However, Wijers et al. (1997) 
also obtain a best-fit GRB standard-candle luminosity of 1.6 x 10 52 erg s~ x (assuming = 1, 
Hq = 50 km s _1 Mpc -1 ), which is ~ 20 times brighter than that obtained for a non-evolving 
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source population. In the following we will consider the best-fit rates for both GRB sources which 
are evolving (following the cosmic star formation history) or non-evolving. 

Next, let us consider a population of transient sources, which are standard candles in peak 
flux and are characterized by a comoving rate per unit volume R{z). The observed number of new 
events per unit time brighter than flux F v at observed wavelength A for such a population is given 
by 

/■ 2 max(F„;A) 

N(F U ;X)= R(z)(l + z)-\dV c /dz)dz, (4) 



o 

where z ma , x (F u , A) is the maximum redshift at which a source will appear brighter than F u at an 
observed wavelength A = c/u, and dV c is the cosmology-dependent comoving volume element. 
The above integrand includes the (1 + z) reduction in the apparent rate due to the cosmic time 
dilation. For Q = 1, A = 0, the comoving volume element covered by solid angle df2 and redshift 
interval dz is given by 



c 



3 (1 + z - y/TTz~) 2 



dv ' = i U) (i + ,)'/» dmz - (5) 

Equation (HJ) is appropriate for a threshold experiment, such as BATSE, which monitors the 
sky continuously and triggers when the detected flux exceeds a certain value, and hence identifies 
the most distant sources only when they are near their peak flux. For search strategies which 
involve taking a series of "snapshots" of a field and looking for variations in the flux of sources in 
successive images, one does not necessarily detect most sources near their peak flux. In this case, 
the total number of events {not per unit time) brighter than F v at observed wavelength A is given 
by 

roc 

N(F V ;X)= R(z%(z;F„,\)(dV c /dz)dz, (6) 



o 

where t+(z; F u , A) is the rest-frame duration over which an event will be brighter than the limiting 
flux F v at redshift z. This is a naive estimate of the so-called "control time"; in practice, the 
effective duration over which an event can be observed is shorter, owing to the image subtraction 
technique, host galaxy magnitudes, and a number of other effects which reduce the detection 
efficiency (Pain et al. 1996). 

In applying equations ([[])-([]), we use the star formation rate p s (z) determined by Madau 
(1997) and convert cosmic time to redshift according to the standard cosmological formulae. For 
the ages of the SNIa progenitors we use the calculations by Schaller et al. (1992). To determine 
the maximum redshift z m3X and the effective duration i* for both SNe and GRBs, we use the 
solar-metallicity, time-dependent, spectral models of Eastman et al. (1994) for SNell, the template 
lightcurves of Riess, Press, & Kirshner (1996) together with the spectral models of Nugent et al. 
(1997) for SNela - extended out to 2.2 /jm with a simple blackbody fit, and the broken power law 
model of Waxman (1997) for GRB afterglows. In all cases, we truncate the emission spectrum at 
rest wavelengths shorter than ~ 10 3 A, to reflect the absorption by HI beyond the Lyman-limit 
in the local environment of the source and the intergalactic medium. For the afterglows, we 
modify Waxman's choice of parameters after peak luminosity so as to match the faster decay slope 
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observed for GRB970228 (Galama et al. 1997; Fruchter et al. 1998) and GRB970508 (Galama et 
al. 1998), namely L u oc t . We also introduce a cutoff time t cut beyond which we truncate the 
afterglow flux so that the source counts will converge (see §3 for the effect of varying t cut ). The 
resulting GRB spectrum as a function of time is 



ay 



erg s x Hz , 



(7) 



o 



-1.1 for t X < tday < t 



''CUt ) 



and 



L v {t) = 2.94 x 10 dU r? 
where t day = (t/day), a = 0.5 for i day < t\ 

t\ = (1.41 x 10 14 Hz/V) 2 / 3 days = (A/2.13/xm) 2 / 3 days. The parameter r\ reflects the 
choice of a typical GRB redshift, with rj = 1 for a non-evolving burst population - corresponding 
to a typical redshift z ~ 1. On the other hand, Wijers et al. (1997) derive a characteristic 7-ray 
luminosity which is ~ 20 times higher than that obtained for a non-evolving population. The 
afterglow luminosity depends on the total fireball energy (Waxman 1997), which is proportional to 
the product of the 7-ray luminosity and the rest-frame duration. Taking all of these dependences 
into account, we derive r/ ~ 10 for a population which traces the star formation history. We express 
the afterglow emission in terms of the luminosity per unit frequency, since detector sensitivities 
are often expressed in Janskys (1 Jy = 10 -23 erg s _1 cm~ 2 Hz -1 ). From equation (0), we have for 
GRB afterglows 



U{z;F v ,X) = min{£ cut , A(z;F u ) °- 9 t x/ ( 1+z) } - A(z; F„) 2 t x/ ^ +z) , 



A(z-F v ) = 



4.0 x W z 7] 
(1 + *) 



2^-1 



El 
Jy 



Gpc 



(8) 



where D^(z) is the cosmology-dependent luminosity distance, D\ J (z) = (2c/Hq){1 + z — \Jl + z) 
for Q = 1 and A = 0. 



3. Results 

We apply equation (Q) to determine the average rate at which GRB afterglows, SNela, and 
SNell should be detected in a field of 1 deg 2 as a function of limiting flux. Figure 1 shows the 
results for different spectral bands, namely the K (A c g = 2.2 /um), R (A e g = 7000 A), B (effective 
wavelength, A efT = 4400 A), and U (A eff = 3650 A) bands. The 5-band rate of ~ 4800 SNell 
yr _1 deg~ 2 brighter than 1 nJy is slightly lower than a recent prediction (Madau 1998) of 25-35 
SNell yr _1 per 4' x 4' field of view of the Next Generation Space Telescope (NGST, Mather & 
Stockman 1996), where the range of values corresponds to a range of assumptions about the 
amount of dust extinction between the observer and the sources. Our slight underestimate of 
~ 21 SNell yr _1 per NGST field results from the fact that we do not correct the star formation 
rate for the effects of dust extinction. Our predictions from equations (||)-(|2|) for the local (z = 0) 
SNIa and SNII rates are consistent with recent observations (Cappellaro et al. 1997; Sadat et al. 



- 6 - 



1998). The SNIa rate we obtain at a redshift z ~ 0.4 agrees to within statistical uncertainties with 
the rate determined from observations of high-redshift SNela (Pain et al. 1996). The Type la 
supernova counts overtake the Type II counts at relatively bright fluxes, since SNela are an order 
of magnitude brighter than SNell and hence dominate the counts in a shallow magnitude-limited 
sample, even though their absolute rate is a fraction (~ 30-50%) of the Type II rate. For both 
types of supernovae, we truncate the spectrum at rest wavelengths shorter than ~ 10 3 A due to 
the expected galactic and intergalactic absorption by neutral hydrogen. 

The predicted GRB afterglow rate in Figure 1 is nearly independent of the observed band, 
since the peak GRB flux in equation (0) is the same at all wavelengths. The sharp "knee" 
at 10 _4 -10™ 3 Jy reflects the edge in the spatial distribution of luminous bursts due to the 
non-Euclidean geometry of the Universe and the sharp decline in the star formation rate at z > 4. 
We extrapolate the star formation rate using a power-law dependence on (1 + z) at z > 4. Our 
extrapolation is consistent with the upper limit on p s (z) around z = 5.5, based on the lack of "V 
dropout" galaxies in the Hubble Deep Field (Madau 1996). Production of bright afterglows, of 
the type observed in GRB970228 or GRB970508, requires an ambient gas density of ~ 1 cm -3 
(Waxman 1997), which is typical of the interstellar medium of disk galaxies. Hence, we expect 
significant Lyman- limit absorption due to HI in the local burst environment, and so we truncate 
the afterglow spectrum at rest wavelengths < 10 3 A. This explains the variation with observed 
wavelength of the rate of faint afterglows. Note that at all observed wavelengths, the afterglow 
rate becomes comparable to the Type II SN rate at magnitudes ~ 18 — 20. This means that if 
afterglows occur with the same frequency as the 7-ray selected GRBs, then they should turn up 
with comparable numbers to supernovae in searches of this depth. 

Figure 2 shows the results of applying equation (||) to obtain the number of events that are 
detectable in a single "snapshot" , taking into account the effective length of time during which 
each type of event should be visible at a given sensitivity. Our prediction for the -R-band counts 
of SNela is consistent with the results of the high-redshift survey of Pain et al. (1996), who found 
three Type la supernovae in the range 21.3 < R < 22.3 across a 1.73 deg 2 area on the sky. Our 
predicted rate is somewhat higher, due to the fact that our naive estimate of the "control time" 
t* assumes 100% efficiency over the duration of the event, while the actual efficiency (and the 
corresponding value of £*) are lower in a real search. Afterglows are less abundant relative to 
supernovae in Figure 2 as compared to Figure 1, because the timescale for the decline of their 
lightcurve around peak flux is an order of magnitude shorter [cf. Eqs. (|j) and (g)]. 

The lower and upper solid lines in Figures 1 and 2 show the predicted afterglow rates with 
and without evolution. The two lines differ substantially at faint fluxes and coincide to within a 
factor of a few at bright fluxes. The coincidence at the bright end results from the fact that the 
adopted GRB luminosities and rates were chosen so as to reproduce the same flux distribution 
of BATSE-selected GRBs. However, at faint fluxes which require sensitivities greater than that 
achieved by the BATSE experiment, the counts converge to different values. The no-evolution 
model predicts many more bursts at redshifts low enough such that the observed frequency falls 
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below the Lyman limit at emission. 

In order to apply equation (||) to the afterglow population, we needed to assume that the 
afterglow emission truncates after some time t C ut, since otherwise the power-law lightcurve implied 
by equation (|7|) would lead to an untruncated power-law in the number counts, and a divergence 
of the number of afterglows at faint fluxes. Observations indicate that GRB afterglows continue 
their power-law decline, L u oc t , over a considerable fraction of a year (Galama et al. 1997), 
and so in Figure 2 we assumed a value of t cut = 1 year. The effect of varying this choice of i cu t 
is illustrated in Figure 3; the position of the number count "knee" shifts to fainter fluxes when 
t cu t is increased since a larger i cut corresponds to a larger effective volume being surveyed. The 
uncertainty concerning t cu t translates to an uncertainty in the expected faint counts. By contrast, 
the afterglow rate predictions displayed in Figure 1 are effectively independent of the detailed 
shape of the lightcurve. 

4. Conclusions 

We have predicted the rate at which GRB afterglows should be detected in supernova searches 
as a function of limiting flux. Requiring that the GRB population would evolve according the 
cosmic star formation history reduces the number counts of afterglows at faint fluxes by ~ 2 
orders of magnitude relative to the no-evolution case (see the solid lines in Figs. 1 and 2). Our 
main result is that the average detection rate of new afterglows (in both evolutionary cases) is 
comparable to that of supernovae for a continuous search with a magnitude limit of K = 18 (cf. 
Fig. 1), assuming that every GRB produces a bright afterglow of the type observed for GRB970228 
or GRB970508. However, the number of detectable afterglows in a "snapshot" search (cf. Fig. 
2) could be an order of magnitude smaller than the number of detectable supernovae, since the 
afterglow emission declines around its peak flux an order of magnitude faster. The actual number 
of afterglows in real surveys might be even lower since only a fraction of all 7-ray selected GRBs 
detected in 1997 were observed to give rise to bright afterglows. (One theoretical interpretation 
of this result could be in terms of the proportionality between the predicted afterglow luminosity 
and the square root of the ambient gas density (e.g., Waxman 1997). The distribution of afterglow 
luminosities is expected to reflect the range of gas densities in GRB environments and might span 
several orders of magnitude.) The reduction factor due to the small appearance probability of 
bright afterglows might be counteracted by a possible enhancement factor due to beaming. If the 
initial gamma-ray emission is beamed to an angle ~ I/70 where 70 is the initial Lorentz factor, 
while the afterglow radiation is produced when the fireball has been decelerated to a modest 
Lorentz factor 7 a , then the number of optically-selected afterglows could be greater by a factor 
~ (7o/7a) 2 than the number of 7-ray selected GRBs (Rhoads 1997). Since typical bursts have 
70 <; 10 2 ~ 3 (e.g., Fenimore, Epstein, & Ho 1993; Woods & Loeb 1995) while the optical afterglow 
emission occurs at 7 a ~ 10 1-2 , this could boost the expected afterglow rates by up to four orders 
of magnitude. Existing supernova searches (e.g., Garnavich et al. 1998; Perlmutter et al. 1998) 
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reach an effective magnitude limit R ;> 22 (Kirshner 1998), and could already place meaningful 
constraints on the beaming of the 7-ray emission. Future searches might also constrain the 
appearance probability of the afterglow emission and the evolution of the GRB population. 

The Next Generation Space Telescope, with its detection sensitivity of ~ 10~ 9 Jy at 1-3.5/um 
(Mather & Stockman 1996), will be able to see GRB afterglows out to redshifts z ;> 10, if they 
exist^J. Afterglows might be the brightest sources at these redshifts aside from early quasars 
(Haiman & Loeb 1997b). Based on the expected supernova rate at z ^> 5 (Miralda-Escude & Rees 
1997), there might be one detectable afterglow per year per ~ deg 2 in the no-evolution model. The 
stretching of the observed lifetime of afterglows from high redshifts and the potential enhancement 
factor for the afterglow rate due to beaming, might bring the detection of afterglows closer to the 
realm of feasibility. The discovery of high redshift afterglows could then be used to trace the star 
formation history out to the "dark ages" of the Universe when the first stars and quasars formed 
(e.g., Rees 1996; Haiman &: Loeb 1997a,b; Loeb 1997). The optical afterglow radiation could also 
be used as a means of identifying the reionization redshift through the spectral location of the 
Gunn-Peterson absorption trough which is produced by the neutral intergalactic medium prior to 
reionization. 

We thank Bob Kirshner for useful discussions and Ron Eastman for providing us with the 
latest model spectra for Type II supernovae. 
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Fig. 1. — Cumulative observed rate N(> F v ) per year per square degree of GRB afterglows (solid 
curves) at four wavelengths, corresponding to the K, R, B, and U bands. The lower solid curve 
assumes the best-fit rate and luminosity for GRB sources which trace the star formation history 
(Wijers et al. 1997), while the upper solid curve assumes the best-fit values for a non-evolving GRB 
population. For comparison, we show the counts for SNe la (dashed curves) and SNe II (dotted 
curves) . 
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Fig. 2. — Cumulative number counts iV(> F u ) per square degree for GRB afterglows, with a cutoff 
time of t cut = 1 yr. Notations are the same as in Figure 1. 
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Fig. 3. — Cumulative B-band counts of GRB afterglows for various values of the cutoff time t cu t- 
From top to bottom, the curves are for t cut = oo, 10 years, 1 year, and 1 month. 



